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Scaling and commensurate-incommensurate crossover for the d = 
2, z = 2 quantum critical point of itinerant antiferromagnets. 
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PACS 71 . 10. Hf - Non-Fermi-liquid ground states, electron phase diagrams and phase transitions 

in model systems 
PACS 71. 10. Fd - Lattice fermion models (Hubbard model, etc.) 
PACS 64 . 60 . A — Specific approaches applied to studies of phase transitions 
PACS 73 . 43 . Nq - Quantum phase transitions 

Abstract. - Quantum critical points exist at zero temperature, yet, experimentally their influ- 
ence seems to extend over a large part of the phase diagram of systems such as heavy-fermion 
compounds and high-temperature superconductors. Theoretically, however, it is generally not 
known over what range of parameters the physics is governed by the quantum critical point. We 
answer this question for the spin-density wave to fermi-liquid quantum critical point in the two- 
dimensional Hubbard model. This problem is in the d = 2, z = 2 universality class. We use the 
Two-Particle Self-Consistent approach, which is accurate from weak to intermediate coupling, and 
whose critical behavior is the same as for the self-consistent-renormalized approach of Moriya. De- 
spite the presence of logarithmic corrections, numerical results demonstrate that quantum critical 
scaling for the static magnetic susceptibility can extend up to very high temperatures but that 
the commensurate to incommensurate crossover leads to deviations to scaling. 



There are strong indications that quantum critical 
points, i.e. critical points at zero temperature, influence 
the physical properties of materials at surprisingly high 
temperature. But the precise region of temperature over 
which this influence is felt is currently not well understood. 
In solvable models of quantum critical behavior, [1] power 
law scaling and universality associated with quantum crit- 
icality were found up to temperatures of order J/2 where 
J is the exchange constant. That is in sharp contrast with 
classical critical points where scaling is usually observed 
only in a very narrow range around the critical point. The 
importance of quantum critical points [2] has thus come 
to the fore in the study of numerous materials, includ- 
ing high-temperature superconductors and heavy-fernrion 
materials where quantum phase transitions and power law 
scaling are observed. [3] 

One particularly relevant case in this context is that of 
itinerant electrons undergoing a paramagnetic Fermi liq- 
uid to spin-density wave (SDW) transition in two dimen- 
sions. The Hubbard model is the simplest microscopic 
model that contains this physics. There is no analog of 



the Ginzburg criterion that allows us to determine the pa- 
rameter range where the influence of the quantum critical 
point is important. In that regime, temperature acts like 
a finite-size cutoff for the correlation length £. In this 
paper, we quantify the range of temperature where quan- 
tum critical scaling is observable in this model, in other 
words we find out whether details of the Fermi surface 
(that lead for example to commensurate-incommensurate 
(C-I) crossovers), logarithmic corrections, or interaction 
effects, lead to sizable deviations from quantum critical 
behavior at finite temperature. 

For this problem, the dynamical critical exponent z 
is equal to two and the corresponding universality class 
(d + z = 4) at the upper critical dimension is ill under- 
stood. [2-4] In particular, the standard Hcrtz-Millis ac- 
tion for quantum critical phenomena is invalid. [5,6] More 
specifically, when the SDW is commensurate at the an- 
tifcrromagnctic wave vector, it has been suggested that 
all the coefficients of the Ginzburg-Landau- Wilson action 
become singular and that the spin susceptibility scaling 
becomes 1/T V with n < 1. [5]. The generic case where the 
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SDW is not commensurate should not have these singu- 
larities. 

An alternative approach is the self-consistent renormal- 
ized theory of Moriya. This theory includes logarithmic 
corrections. [3, 7] However, it is not adequate to make 
quantitative predictions for deviations from quantum crit- 
ical effects in the Hubbard model since it necessitates phe- 
nomenological constants as input. In addition, it does not 
satisfy the Pauli principle. In a theory that satisfies the 
Pauli principle, there is an interaction independent sum 
rule on spin and charge susceptibilities [8] that should be 
enforced and, in addition, the local moment, (<Sf), with 
S z the z component of the local spin, cannot exceed h 2 n/A 
when the filling n satisfies n < 1 and fi 2 (2 — n)/4 when 
n > 1. There is nothing that imposes these constraints in 
the theory of Moriya. 

Method and model: We use the non-perturbative Two- 
Particle Self-Consistent (TPSC) approach [8]. This ap- 
proach respects the Pauli principle, the Mermin- Wagner 
theorem and conversation laws. It also contains quantum 
fluctuations in crossed channels that lead to Kanamori- 
Briickner screening. [9] It is valid in the weak to inter- 
mediate coupling regime (U < 6i) and not too deep in 
the renormalized classical regime where a pseudogap is 
observed. Numerical results obtained from TPSC in its 
domain of validity are extremely close to the numerically 
exact solution obtained (barring statistical errors) with 
benchmark Quantum Monte Carlo calculations on the 
Hubbard model. [8-13]. The approach gives a satisfactory 
description of the pseudogap in electron-doped cuprates 
in a wide doping range. [14, 15] It has been shown to be 
in the N = oo universality class, where N is the number 
of components in the O(N) vector model. [16] Since we 
are looking for deviations from universality and the the- 
ory has been benchmarked in non- universal regimes, we 
argue that our results are reliable for this question, even 
though we cannot claim to be completely accurate in the 
N = 3 regime. Nevertheless, we will demonstrate that 
TPSC has the same critical behavior as Moriya theory 
and hence has the same logarithmic corrections. These 
logarithms have the same functional form as those of the 
renormalization group asymptotically close to the quan- 
tum critical point, but in TPSC and in Moriya theory the 
mode-mode coupling term does not flow, hence the cor- 
rections may differ in the details from the renormalization 
group. [3] Quantum critical behavior of the susceptibility 
and of the self-energy in the closely related spin-fermion 
model has been discussed by Abanov et al. [17] 

We study the t—t' — U two-dimensional Hubbard model 
on the square lattice at weak to intermediate coupling, 

H = - Y^ *»,.? (4,<r c J> + h - c -) + U ^2 n ht n i,l C 1 ) 

where tij are the hopping integrals, i,j are the site index, 
a is the spin label, c\ a and Ci y(J are the particle creation an 
annihilation operators. Each doubly occupied site costs an 



energy U and rii ll7 = c\ a Ci. a . The units are such that h = 
1, ks = 1 and lattice spacing is unity. All the numerical 
results are presented in units where t = 1. The dispersion 
relation is written as: 

£k = — 2t(cos{k x ) + cos(ky)) — At' cos(k x ) cos(k y ). (2) 

We concentrate on the behavior of the spin susceptibil- 
ity In TPSC, the retarded spin susceptibility x(q, u>) is 
written as: 

Xo(q,w) 



x(q,w) 



i 



(3) 



-J £ xo(q s w) 

where XoiQ.^) is the retarded Lindhard function at wave 
vector q and angular frequency to. The effective spin in- 
teraction U sp is evaluated without adjustable parameter 
using the ansatz [8, 9] 

U(n 7 ni) = U sp (n^(ni) (n < 1) (4) 

U((l-n 1 )(l-n l )) = U sp ((l-n 1 ))((l-n l )) (n > 1) (5) 

with the local-moment sum rule that follows from the 
fluctuation-dissipation theorem 

dw f°° d 2 q 2 

2?r 



n-2(n T nj.) 



(2nY 1 



e -u>/T 



X"(q,w) 



(6) 

where x"(q, w) = Imx(q, u>), T is the temperature, and 
(n-j-nj.) double occupancy. We dropped the site index using 
translational invariance and we used the Pauli principle to 
write 

S 2 ^((n^-n i ) 2 )=n-2(n m ). (7) 

All the numerical results below are obtained using the 
Matsubara frequency version of equations (2) to (6) with- 
out any approximation, hence they are valid at arbitrary 
distance from the quantum critical point. Before proceed- 
ing, we show however that the quantum critical behavior 
of TPSC is the same as that of the self-consistent renor- 
malized theory of Moriya and we discuss conditions for 
scaling. 

Analytical results near the quantum critical point. 
When the correlation length is large, one can expand the 
denominator of the TPSC spin susceptibility around the 
wave vectors q^ where the maxima in xo occur to obtain 



x"(q,w) 



lo/Tq 



Uspt,Q 



2 /_^, 



(e- 2 + (q-q«) 2 ) 2 + (c/r ) 2 ' 



(8) 



Defining U m f = 2/xo (qi,0) as the value of the interaction 
at the mean-field SDW transition, the other quantities in 
the previous expression are 
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5U 
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u, 



4o1 su 



Urnf Usp? 



i 



d 2 Xo (q,0) 



2xo (0,qi) dq 2 
1 dxo(cii,io) 



£oXo(qi,0) d(iu>) 



(9) 

(10) 

(11) 

(12) 
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In the expression for the spin susceptibility, the denomi- 
nators are expanded around each of the four incommensu- 
rate wave vectors, or only around the (tt, tt) wave vector 
depending on the situation. We checked explicitly that 
higher powers of (q — q^) do not improve the description 
of the C-I crossover and are not relevant. 

To determine the quantum critical behavior, one sub- 
tracts the self-consistency condition Eq.(6) for a value of 
temperature and filling close to the quantum critical point 
from the same equation evaluated at that critical point 



s A - si 



dui 

TT 



<Pq 

(27T) 2 



j/T _l 



-x"(q,w) 



+x"(q,«)-x£(q,«)] 



(13) 



In the above expression x" (li w ) i s evaluated at the quan- 
tum critical point where £ ~ 2 = 0. (From now on, a sub- 
script c means that the quantity is evaluated at the quan- 
tum critical point.) One then performs the integrals over 
momentum in a circular domain with cutoff qs and then 
the frequency integrals. To write the final answer, it is 
useful to follow Moriya et al. [18] and to define 



Tn = 



r g| 

2tt ' 



and dimcnsionless measures of £ and T 



r 2 

v = ~r ; T 

1b 



T 



(14) 



(15) 



The definition of Tq Eq.(12) and the fact that £o and q^ 1 
are both of the order of the lattice spacing shows that To is 
a temperature of the order of the Fermi energy. With these 
definitions and a single maximum in the susceptibility, the 
self-consistency expression takes the form 



2/(1 -In y) = Vo + - 

7T 



y_ + i_ 

T T 



(16) 



where terms of order y 2 have been neglected on the left- 
hand side. We defined 



Vo 



UspZ,Q 
Tn 



{s 2 - s 2 ) 



(17) 



and obtained <\> (x) from the second Binet log gamma for- 
mula [19] 



cf>(x) 



1 ( z 

2 / dz^ arctan — 

e 2vz - 1 V x 



o 



(18) 



lnT(x)- ( x- -) \nx + x- -ln(27r)(19) 



with i?e[z] > and T (x) Euler's gamma function. The 
quantity yo in Eq.(17) measures the deviation from the 
quantum critical point. Apart from the logarithm, the 
self-consistency relation Eq.(16) has the same functional 



form as Eq.(2.8) in Ref. [18] Logarithmic corrections for 
that theory are mentioned without proof in Ref. [7] . 

For large local moment, S 2 > S 2 , there is an SDW 
ordered ground state and yo < 0. The case yo > cor- 
responds to the Fermi liquid ground state and yo = to 
the quantum critical point. The full filling and tempera- 
ture dependence of yo is found from the definitions of U sp 
and £ -2 . For example in the hole doped case, defining 
An = n — n r . we have 



yo 



UmfZl 



1 



T„ 






+ 1 



An 



U, 



mf 



u. 



2U 






1 



c,mf 2 



2U 



(20) 

Thus, yo depends on y = £ 2 / ' q\ but in the critical regime 
y C 1 and yo -C 1 so we can neglect terms of order yyo. 
The quantity yo can then be written in the form 



yo 



(aAn + bT) 



(21) 



where a is a positive number. 

The various limiting solutions for the critical behavior 
of the dimensionless correlation length can be obtained 
from the self-consistency condition Eq.(16) as follows. [18] 
For yo < 0, one must take the limit y — > first, then 
<j> (^) — 4> (^ + i) ~ — i In (-^r-) and since y is exponen- 
tially small, y — y In y can be neglected on the left-hand side 
leading to y ~ (r cxp(27ryo/r))/27r. This is the renormal- 
ized classical regime where the correlation length grows 
exponentially At the quantum critical point yo = 0, the 
same limit of the <j) functions applies and one must find 
the solution of — ylny ~ — ^-lnf-^J which is approxi- 
mativcly y ~ rln (|1iit|) / |hir|, as in the renormalization 
group. [7] Finally, in the Fermi liquid regime, yo > 0, 
the correlation length (and hence y) is finite so the r — > 
limit must be taken first and <f> (£) - <f> (^ + i) ~ r/ (12y) 
which yields y ~ yo + O (r 2 ) . At r = on the Fermi liquid 
side, there are logarithmic corrections to the dependence 
of y on yo asymptotically close to the quantum critical 
point since — ylny ~ yo whose approximate solution is 

y— -yo/hiyo. 

In all regimes where lny in the self-consistency Eq.(16) 
can be neglected (large T) or replaced by a constant in 
the temperature range of interest, one can write 



y =F /An l 



where the scaling function F is the solution of 



cF = 



.'/o 



1 

7T 


<P{F)- 


-*(, + i)] 



(22) 



(23) 



with c = 1 — In y t , y t being the typical value of y in the 
range of temperature under study. We have already dis- 
cussed limiting cases of F above. We demonstrate numer- 
ically below that in the range O.Olt < T < t logarithmic 
corrections are negligible so that scaling holds to an excel- 
lent approximation, except at the C-I crossover. 
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Scaling function. When the explored temperature range 
is limited on a logarithmic scale, or when T is large, log- 
arithmic corrections can be neglected. In addition, in the 
limit where r is much smaller than y/r, the scaling func- 
tion F in Eq.(22) depends only on An/r since we are in the 
limiting case <p [F + -) — > (ft (oo) = in the equation that 
defines F, Eq.(23). This case occurs when the ground state 
is paramagnetic, tjq > 0, or above the crossover line to the 
renormalizcd classical regime that occurs when yo < 0. In 
such cases, near anyone of the maxima located at q^, the 
quantity y = £~ 2 / ' q 2 B scales as tF(Au/t, oo) so the spin 
susceptibility Eq.(8) as a function of an arbitrary scale 
factor s obeys the scaling relation 

x(T,An,|q-<b|,w)-i2 

Xi(8 1 ' v T,8*'"An,8\< l -q i \,a z w) (24) 



<rtl v 



where the exponents have values 7 = 1, v = 1/2, z = 2 
and (/) = 1. In the above equation, R will not be important 
only if the incommensurate peaks are much narrower in 
momentum space than the inverse correlation length. Let 
lo = for now and drop the dependence on that variable. 
Following the above discussions on the behavior of the 
correlation length, the susceptibility xi on the right-hand 
side of the last equation should be, within log corrections, 
a universal function of its arguments but with the overall 
scale of each argument and of xi non-universal. Setting, 
q = |q— q; | = 0, cj = and choosing s such that Ts x l v = 1 
we find, 



1 / An 
x (T,An,Q,0) = -X( — 



R. 



(25) 



where the scale of the function X defined by this equation 
and an overall prefactor in front of the argument are not 
universal. X is the quantity we will focus on, but we note 
in passing that the general form Eq.(24) with the given 
exponents implies lo/T scaling for the q integrated suscep- 
tibility. [20] Non- universal factors such as U sp , £0 and r 
that enter the spin susceptibility can have some tempera- 
ture and filling dependence in TPSC that can in principle 
lead to deviations to scaling. In renormalization group 
language, these dependencies are the irrelevant variables 
whose importance we are trying to gauge to delimit the 
scaling regime. 

Commensurate-incommensurate crossover. In a strict 
sense, the value of q should be fixed at q = q i (T = 0) 
to check quantum critical scaling. However, q^ itself de- 
pends on temperature in general. At high temperature q^ 
equals Q = (tt, it), becoming incommensurate at low tem- 
perature. The susceptibility there shows four symmetry 
related peaks for the model we consider. [21] The value 
of q, (T) clearly depends on details of the Fermi surface 
and is thus non-universal. The above scaling form Eq.(24) 
nevertheless suggests that scaling in the T, An plane as 
in Eq.(25) should occur when q = q i (T) . It is not how- 
ever possible to define q^ (T) in the C-I crossover regime. 




iog 10 (T) 



Fig. 1: (Color online) log 10 (x) as a function of log 10 (T) for 
U = 4t and t' — (open circles) at n c = 0.84 and for the 
corresponding non-interacting susceptibility at U — (open 
squares). The vertical dashed lines indicate the commensurate 
to incommensurate crossover region. The red dash-dotted line 
has a slope — 1. The black vertical line simply indicates T — t. 



In that regime, incommensurate peaks necessarily over- 
lap since the second derivative of x vanishes at q = Q 
when the crossover begins, reflecting the fact that there 
is a broad maximum at Q that is splitting into four over- 
lapping peaks. R in the general scaling function Eq.(25) 
is not negligible in the C-I crossover region. On general 
grounds then, we expect deviations to scaling there. One 
may think that a better strategy to prove scaling is to mea- 
sure the correlation length £ as a function of T an An, but 
£ cannot be determined in the C-I crossover regime for the 
same above reasons. 

From now on, we thus look for scaling with the suscepti- 
bility evaluated at its maximum, x(T, An, |q m ax — qi|, 0). 
This is a well defined quantity experimentally and far from 
the C-I crossover we will have q max = qi (T) . 

Numerical results: Let us first verify the scaling at the 
quantum critical point An = 0. Fig. (1) shows a log- 
log plot of both the interacting (open circles) and nonin- 
teracting (open squares) susceptibilities as a function of 
temperature for two different sets of parameters. 

For temperatures larger than hopping t, one obtains 
trivial \/T scaling for both the interacting and non- 
interacting susceptibilities. While the non-interacting sus- 
ceptibility flattens at lower temperature, the interacting 
susceptibility shows quantum critical 1/T scaling down to 
the lowest temperature we could reach, namely T = O.Olt. 
We will see that the 1/T scaling at T > 1 that comes from 
the non-interacting susceptibility does not obey the scal- 
ing equation Eq.(25). It is also clear from Fig. (1) that 
deviations to scaling occur in the C-I crossover regime de- 
limited by the vertical red lines. It is remarkable however 
that the same straight line fits both the commensurate 
and the incommensurate regimes. This suggests that non- 
universal scale factors are very similar on either sides of 
the commensurate-incommensurate transition. The slight 
upward curvature at the lowest temperatures is not incon- 
sistent with effects of logarithmic corrections. 
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To verify the full scaling Eq.(25), we plot Tx as a func- 
tion of \An\/T on a log- log plot in Fig. 2. We take values 
of n on the Fermi liquid side of n c . For a given band 
structure and interaction, it is only when one has found 
the correct values of the critical n = n c that all the curves 
for different fillings and temperature collapse on the same 
curve. We found, when t' = 0, that n c = 0.926, 0.840 and 
0.795 for U = 2t, At and 6i respectively and n c = 1.180 for 
the electron-doped case with U = 6t, t' = — 0.054. More 
values can be found in the thesis which is the basis for all 
the results of the present paper. [20] 

The straight line of slope —1 at large An/T in Figs. 
2a,c,d. corresponds to the Fermi liquid regime where both 
the susceptibility and the correlation length are tempera- 
ture independent, but diverge as one approaches the quan- 
tum critical point. In that regime, x scales as £ 2 ~ 1/ An 
when logarithmic corrections are negligible. The 1/T 
scaling of x corresponds to plateaus on the left of Figs. 
2a,c,d. The deviations from a plateau come from the C-I 
crossover. To show that the scaling is non-trivial, in Fig. 
2b we do not multiply the susceptibility by T on the ver- 
tical axis. The lined up circles that can be caught by the 
eye correspond to different temperatures for a given filling 
n, the fillings closest to n c being to the left. 

Scale factors depending on band structure and interac- 
tion strength should not influence the shape of the scaling 
function. A simple translation in the T — An plane of the 
curves for different parameters should allow all of them to 
collapse. In Fig. 2d, we show scaling functions for various 
parameters but without translation for non-universal fac- 
tors. One sees that if there were no deviations to scaling 
associated with the C-I crossover in the plateau region, 
simple translation would make all the curves nearly col- 
lapse. This also shows that logarithms do not have a large 
influence on scaling in this temperature range. 

In Fig. 3, data analogous to those in Fig. 2 are repre- 
sented by black open circles and are filtered out near the 
C-I crossover. The missing data is particularly clear in 
Fig. 3b where we do not scale the vertical axis. If Ti is the 
temperature where the crossover occurs for a given doping, 
the data were filtered in the range T { - AT < T < T t + AT 
(AT ~ 0.2i) for densities n c — n < 0.04. For larger values 
of n c — n, the data is sufficiently far from the C-I crossover 
that no filtering is required. The remaining data are those 
beyond the C-I crossover both above (commensurate) and 
below (incommensurate) Tj. One sees that a plateau is 
recovered (black open circles) for all three values of the 
interaction strength appearing in Figs. 3a,c,d as expected 
in the quantum critical regime. 

We now turn to the high temperature limit of the quan- 
tum critical scaling. While the black open circles in Fig. 
3 are for T < t, those for t < T < lOt are represented 
by red crosses. The deviations to scaling for t < T < lOt 
are obvious. Even though the non-interacting susceptibil- 
ity scales as 1/T for T > t as we saw in Fig. 1, it does 
not pollute the scaling associated purely with the quan- 
tum critical point. The latter occurs for T < t, with the 
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Fig. 2: (Color online) a) log 10 (Tx) as a function of log 10 ( -^p- 1 ) 
for U = At and t! — 0. b) Shows the data of figure a) with the 
susceptibility unsealed log 10 (x) c) log 10 (Tx) as a function of 
log 10 ( h_£l J for U — 6t and t' = — 0.05£. d) compares the 

scaled data log 10 (Tx) as a function of log 10 (-^?) for U — it 
and t' — (green crosses), U — &t and t' — (red open squares) 
and U — 6t and t' — — 0.05£ (black open circles). On figures a) 
and c) we can see deviations from the scaling at log 10 h=M < — 1 
that are due to the commensurate incommensurate crossover 
(see text). In these figures, for t ' = we considered hole doping 
with n c — n > 0, the smallest n being n = 0.7. For t' = — 0.05£, 
we took electron doping with the largest n equal to n = 1.3 
and n — n c as argument of the log. 



caveat concerning the C-I crossover. The maximum T for 
scaling, T ~ t, is an important result that applies in the 
weak to intermediate coupling regime we have considered 
here. Clearly quantum critical behavior must disappear at 
U = 0, so there should be some U dependence to the up- 
per temperature cutoff. At the intermediate coupling val- 
ues that we considered, the temperature range over which 
quantum critical scaling is observed should be compared 
to what would have been naively estimated by substitut- 
ing U = At and U = 6t in J = At 2 /U [1], obtaining re- 
spectively J/2 = t/2 and J/2 = t/3. Basically, the upper 
limit of T ~ t is essentially the degeneracy temperature 
for Fermi-Dirac statistics, which is of the same order as To. 
The irrelevant temperature dependencies of all quantities 
are thus on this scale. 

Conclusion. The quantum critical behavior of TPSC 
for the d = 2, z = 2 universality class is the same as 
that of the self-consistent renormalized theory of Moriya, 
hence it includes logarithmic corrections. In TPSC there 
is no adjustable parameter. By explicit numerical calcu- 
lations away from the renormalized classical regime of the 
d = 2 Hubbard model in the weak to intermediate cou- 
pling, we have been able to show that logarithmic correc- 
tions are not really apparent in the range of temperature 
0.01 < T < t and that the maximum static spin suscep- 
tibility in the (T, n) plane obeys quantum critical seal- 
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Fig. 3: (Color online) The black circles are the data for T < t 
and the red crosses for T > t. a) Scaled data log 10 {T\) for U — 
At. b) Unsealed values log 10 (x) for the same U. c) log 10 {T\) 
for U = 6t d) for U = 2t. All four panels are at t' = 0. The 
high temperature limit of the quantum critical scaling is of the 
order T ~ t. 



ing. However, near the commensurate-incommensurate 
crossover, one finds obvious non-universal temperature 
and filling dependence. Everywhere else, the (T, n) depen- 
dence of the non-universal scale factors is relatively weak. 
Strong deviations from scaling occur at temperatures of 
order t, the degeneracy temperature, reflecting the fact 
that the temperature dependence of most irrelevant terms 
is on the scale of the Fermi energy. That high temperature 
limit should be contrasted with J/2 found in the strong 
coupling case. [1] In generic cases the upper limit T ~ t 
is well above room temperature. In experiment however, 
the non-universality due to the C-I crossover may make 
the identification of quantum critical scaling difficult. And 
since the (T, n) dependence of q; is non universal, one may 
encounter cases where this is in practice impossible. 

Electron-doped high-temperature superconductors ap- 
pear as an ideal system to check quantum critical scaling 
since they seem well described by the d = 2 one-band 
Hubbard model at weak to intermediate coupling. [14, 15] 
And experiments [22,23] strongly suggest the presence of 
a quantum critical point in these materials. In the case of 
heavy fermions there are examples of SDW-Fermi liquid 
quantum critical behavior. [3] However, these are multi- 
band systems where there are additional energy scales, 
such as the Kondo coherence scale, so our results would 
apply only in regimes where an effective one-band Hub- 
bard model applies. 
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